Multiflavor and multiband observations of neutrinos from core collapse supernovae 
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It has been proposed that the gamma ray burst - supernova connection may manifest itself in a 
significant fraction of core collapse supernovae possessing mildly relativistic jets with wide opening 
angles that do not break out of the stellar envelope. Neutrinos would provide proof of the existence 
of these jets. In the present letter we calculate the event rate of >100 GeV neutrino-induced cascades 
in km^ detectors. We also calculate the event rate for >10 GeV neutrinos of all flavors with the 
DeepCore low energy extension of IceCube. The added event rate significantly improves the ability 
of km'^ detectors to search for these gamma-ray dark bursts. For a core collapse supernova at 10 Mpc 
we find ~4 events expected in DeepCore and ~6 neutrino-induced cascades in IceCube/KMSNet. 
Observations at >10 GeV are mostly sensitive to the pion component of the neutrino production in 
the choked jet, while the >100 GeV depends on the kaon component. Finally we discuss extensions 
of the on-going optical follow-up programs by IceCube and Antares to include neutrinos of all flavors 
at >10 GeV and neutrino-induced cascades at >100 GeV energies. 

PACS numbers: 96.40.Tv, 97.60. Bw, 98.70.Sa 



I. INTRODUCTION 



Long duration gamma-ray bursts (GRBs) and core 
collapse supernovae are known to be correlated (T], 
MeV photons that give rise to GRBs are thought to 
be produced by accelerated electrons in internal shock 
of jets emitted by the GRB progenitor. Many details 
about GRB phenomenology remains uncertain. These 
GRB jets have very narrow opening angles and very 
large Lorentz boost factors (F > 300.) GRBs are one 
of the plausible candidates sources for the highest en- 
ergy cosmic rays [1, 01 and as such would observable in 
~ 100 TeV neutrinos by km"^ detectors such as IceCube 
and KMSNet 0. 

Only a very small fraction of core collapse supernovae 
result in GRBs. It has been speculated that the rate of 
production of jets in core collapse supernovae may be sig- 
nificantly higher than the rate of GRBs, but that often 
these jets are choked within the supernovae. Evidence 
for this hypothesis is the observed asymmetry in the ex- 
plosion of core collapse supernovae 0, . Very recently 
evidence for mildly relativistic jets (F ~ 1) that did break 
from the progenitor has been presented for type Ic super- 
novae 2007gr and 2009bb Supernovae with hidden 
jets, supernovae with mildly relativistic jets that manage 
to break out and long GRBs could form a continuum class 
of astronomical objects. Neutrinos would provide us with 
information about these hidden jets and early neutrino 
observations could also lead the way in finding and study- 
ing supernovae with mildly relativistic jets that do break 
out. A model has been proposed by Razzaque, Meszaros 
and Waxman (henceforth RMW) [ll|. This model was 



extended by Ando and Beacom (henceforth AB) to in- 
clude kaon production [l^. The RMW/AB spectrum is 
very soft leading to >100 GeV neutrino observations in 
IceCube/KMSNet. A model for neutrino production in 
reverse shocks for both choked and successful relativistic 
jets associated with SN type lb has also been proposed 
nj. While we did not consider this latter model in the 
present letter, the broad conclusions still apply. 

A promising way to search for -^100 GeV neutrinos 
from core collapse supernovae is the optical follow-up 
(l^ . Neutrino multiplets in directional and time coin- 
cidence trigger observations by fast robotic telescopes. 
This method has the advantage of being able to dra- 
matically reduce the intrinsic atmospheric neutrino back- 
ground because the signal search window is only O(lOOs). 
A neutrino multiplet may not be significant on its own 
because the accidental rate due to atmospheric neutrino 
background is O(10/year), but the subsequent observa- 
tion of the rising light-curve of a supernova in spatial and 
temporal coincidence with a neutrino multiplet would be 
highly significant. Optical follow-up programs are al- 
ready in operation by IceCube [l5j and Antares [l6j . 

The objective of this letter is to calculate the event rate 
from core collapse supernovae as observed by ~10 GcV 
detectors like DeepCore [17|. We propose the optical 
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follow-up programs to be enhanced so as to require at 
least one ^100 GeV (minimum requirement to achieve 
good sky localization) and several >10 GeV events. We 
also calculate the event rate on > 100 GeV neutrino- 
induced cascades. A most interesting case is that of 
KMSNet, which promises good angular reconstruction of 
cascades. We propose the optical follow-up programs to 
be enhanced so as to require at least one ~100 GeV 
and one ^100 GcV neutrino-induced cascade. The en- 
hanced sensitivity of these two new modes of observation 
will allow IceCube/KMSNet to provide a more thorough 
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test of the RMW/AB model. 



III. EFFECTIVE AREAS AND EVENT RATES 
IN ICECUBE/KM3NET/DEEPCORE 



II. NEUTRINO PRODUCTION IN CHOKED 
JETS 



The RMW/AB model supposes a choked jet with 
bulk Lorentz boost factor F;, ~ 3 and an opening an- 
gle 9j ~ 0.3. The kinetic energy of the jet is set to 
Ej = 3 X 10^^ erg in analogy to GRBs. Also simi- 
lar to GRBs the time variability of the engine is set to 

~ 0.1 s. Neutrinos are produced in p-p interactions 
via both pions and kaons. The accelerated proton spec- 
trum is assumed to be dNp/dE = E~^. The density 
and energies involved are similar to those of atmospheric 
neutrinos in which muons do not decay, but instead are 
subject to radiative cooling. Thus in the case of neu- 
trinos generated by pion decay, the neutrino flavor flux 
ratio (j)u^'-4'vmn'-4'v.r is 0:1:0. In the case of kaons there 
is a small flux of v^. due to decay, but the AB paper 
does not include neutral kaon contribution. We therefore 
also assume 0:1:0 as the flavor flux ratio. Taking into ac- 
count vacuum oscillations the expected flavor flux ratio at 
Earth is 0.2:0.4:0.4 for both pions and kaons. Note that 
vacuum oscillations were not taken into account in pre- 
vious calculations [ill, Ell ' include them here as they 
are critical to describe Ve and Vt fluxes at Earth. All 
our numbers include vacuum oscillations except where 
explicitly noted. 

The mesons product of p-p interactions have the same 
initial spectrum as the protons, but they are subject to 
hadronic and radiative cooling. This results in two break 
energies above which the meson spectrum is steeper. 
Neutrinos follow the energy spectrum of their parent 
mesons. The neutrino flux resulting from pion and kaon 
contributions can be described as a doubly broken power 
law: 
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A is the flux normalization (for all flavors com- 
bined) such that for pions(kaons) dN^/dE is 5 x 
10-2 GeV-i.cm-2 (5 x 10"^ GeV^cm-^) at e'^^K The 

energies E^J-^ and e'^J-'^ are the energies above which 
hadronic and radiative cooling become relevant respec- 
tively. For pions (kaons) these energies are: e'i^^ = 

30 GeV (200 GeV) and E'i^^ = 100 GeV (20 TeV). The 
values above correspond to a supernova at 10 Mpc with 
default choices of parameters. See RMW[TT| and AB[T2| 
for details. 



To calculate the expected number of events Nobs given 
a neutrino flux dN^/dE the neutrino effective area A^ff 
of a detector is needed: 



obs 



(2) 



The neutrino effective area of IceCube for muon neutri- 
nos averaged over the northern hemisphere and assuming 
equal v and v fluxes has been published by the IceCube 
collaboration [l8| and it is reproduced in FiglT] The neu- 
trino effective area of a detector can be calculated from 
first principles and the procedure for doing so has been 
reported in many publications. In particular (l9j and 
refs. therein have a careful discussion of this subject. 

The calculation of neutrino effective area requires 
knowledge of the muon effective area, which in turn de- 
pends on detailed knowledge of the operation of the de- 
tector. In the literature it is common to assume that 
IceCube has an effective area of 1 km^ for muons above 
a threshold of 100 GeV. This approximation is a good 
one for the hard spectra, such as that expected from 
the sources of cosmic rays. In contrast this approxima- 
tion is a bad one for soft fluxes, such as those from the 
RMW/AB model. For RMW/AB most of the events 
observed have energies that are close to the muon detec- 
tion threshold and thus detailed knowledge of the detec- 
tor functioning (and hence the muon effective area) are 
needed. As similar argument also applies to DeepCore 
and cascade events. 

We have used ANIS [l^] to simulate neutrinos in the 
vicinity of IceCube and DeepCore. We assume that Ice- 
Cube is a cylinder of 1000 m of height and 564 m in 
radius. The center of this cylinder is placed 1950 m be- 
low the ice surface. DeepCore is simulated as a cylinder 
of 350 m of height and 125 m of radius at 2275 m of 
depth below the ice surface. The ANIS simulation takes 
into account the ice/rock boundary. Earth neutrino ab- 
sorption (though this is a very small effect for the soft 
RMW/AB spectrum), neutral current regeneration, etc. 
Our simulations include neutrinos of all flavors. We have 
assumed flavor flux ratio as appropriate for an the in- 
jected fluxes by either pions or kaons and we have taken 
into account neutrino oscillations. We have also assumed 
equal ratios of neutrino to anti-neutrino for a given fla- 
vor. Muon propagation is calculated with MMC [2l| . We 
assume that events are detected for the following condi- 
tions: a) through-going muons or entering muons must 
have a minimum energy when entering the detector cylin- 
der b) the total visible energy (cascades, muons, etc) of 
a contained event must be greater than the same thresh- 
old as through-going events. For IceCube we have chose 
100 GeV as the threshold and for DeepCore 10 GeV. 
IceCube results also apply to KM3Net. We show our 
calculations for neutrino effective areas for IceCube and 
DeepCore in Fig. [T] 
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FIG. 1: The three panels show the neutrino effective areas for i/e (left), (center) and i/t (right). The effective areas shown are 
for IceCube (solid histogram - black online version) and DeepCore (short dashed histogram - blue online version) as calculated 
in this letter. Also shown (long dashed curve - red online version) is the effective area for published by the IceCube 
collaboration. The effective area has been averaged over the northern hemisphere {S: 0° to 90°). The effective area is also 
averaged over neutrinos and anti-neutrinos. The difference between IceCube's effective area for nUmU as calculated in this letter 
and as published by the IceCube collaboration is due to threshold effects. In this letter's calculation we have simplistically 
assumed that all muons with i 100 GeV energy that go through IceCube are detected. For details see section Hill 



Note good agreement between our calculations and Ice- 
Cube published effective for >TeV. As described above, 
for low energies, detailed simulation of the detector per- 
formance becomes important. For DeepCore our simu- 
lation is in good agreement with the results published 
by the IceCube collaboration The IceCube collab- 
oration has not published an effective area for neutrino- 
induced cascades. Therefore our results concerning cas- 
cades will have the largest uncertainty. 

Performing the integral on Eq. [5] using IceCube's pub- 
lished neutrino effective area with the RMW/AB flux for 
a core collapse supernova at 10 Mpc confirms previous es- 
timates of expected i^^ events Using IceCube' pub- 
lished I'fi effective area we obtain an expectation of 11.2 
events (28 if we don't assume oscillations as AB did.) 
The comparison remains equally valid at high energies 
for which threshold effects are less relevant. Using our 
calculation of the effective area we obtain 3.3 events. 
Our discrepancy with AB and the IceCube published ef- 
fective area can be attributed exclusively to threshold 
effects near 100 GeV. 

The all-flavor expectation for DeepCore is 4 events. 
The neutrino-induced cascade signal expected in IceCube 
is 6.1 events. DeepCore events are broken down by flavor 
to 0.4 due to Ve, 3 are due to u^^ and 0.6 due to Vr- 
Cascade events break down to 2 due to Ve, 1-3 due to 
(via neutral current interactions) and 2.8 due to Vr- 

Figure shows the observed event spectra for all three 
neutrino flavors in DeepCore due to a reference super- 
nova at 10 Mpc that follows the RMW/AB model. Fig- 
ure [3] is similar but neutrino-induced cascades in a km"^ 
detector. It is clear that the pion component makes the 
most significant contribution to the observations in Deep- 
Core. In contrast events expected in IceCube, both 
and cascades, are mostly due to kaons with the exception 
of energies near to the detector threshold. 



IV. ATMOSPHERIC NEUTRINO 
BACKGROUND 



In DeepCore 10^ atmospheric neutrinos of all flavors 
over 27r sr are expected per year [l^. IceCube has an 
atmospheric muon neutrino rate of 10^ events per year 
over 27r sr 0]. IceCube has a rate of 2 x 10^ atmospheric 
neutrino- induced cascades per year over 47r sr[22l|. 

When searching for neutrino transients, background is 
significantly reduced by using a narrow time window and 
by pointing in a specific direction. Km^ neutrino detec- 
tors have angular resolutions of 0(1°) H^]. Due to the 
optical properties of ice, IceCube is not expected to be 
able to reconstruct the direction of cascades well, but it 
is anticipated that KM3Net will be able to reconstruct 
cascades with 5° resolution [2^. It is also expected that 
DeepCore will reconstruct the direction of events. For 
charged current j/^ events above ^10-20 GeV DeepCore 
will have a resolution of 10-15°. For cascade- like events in 
DeepCore the angular resolution is expected to be ~ 30° 
but the minimum energy at which this is feasible is still 
unknown. DeepCore will also have the ability to separate 
track events (from C.C. interactions) from cascade-like 
events for energies greater than 10 GeV [l^]. In our cal- 
culations below we assume, perhaps pessimistically, an 
angular resolution of ~ 30° for all flavors in DeepCore. 
The appropriate time search window is set by the size 
of the star: Rt,/c\ight ~ 100s). In fact the radii of pro- 
genitors just prior to core collapse are known with poor 
detail. Candidate progenitor stars have been identified 
for about a dozen of core collapse supernovae, including 
SN1987A (type IIP, i?* ^ 10" cm) ^| and SN2008D 
(type lb, i?* - 10" cm) But SN1987A is unusual 

and type II progenitors are more frequently red giants 
with i?* as large as ~ 10^"^ cm. 

The rate of accidental atmospheric neutrino multiplets 
with N events of channel a (e.g. >100 GeV i^^) and M 
(e.g DeepCore) events of channel b is: 




FIG. 2: The panels show the expected signal above a given energy for DeepCore for Ve (left), (center) and i^t (right). The 
solid (black in the online version) histogram is the total expectation of the RMW/AB model, the long dashed histogram (blue in 
the online version) is the pion contribution and the short dashed (red in the online version) histogram is the kaon contribution. 




FIG. 3: The panels show the expected signal above a given energy for cascades in IceCube/KMSNet for Ue (left), i/^i, (center) 
and v-r (right). The solid line (black in the online version) is the total expectation of the RMW/AB model, the long dashed line 
(blue in the online version) is the pion contribution and the short dashed line (red in the online version) is the kaon contribution. 



V. DISCUSSION 
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where Ra and Rh are the atmospheric neutrino back- 
ground rates for channels a and b, Q.a and f2f, are the an- 
gular areas search areas and AT is the temporal search 
window. Note that eq. |3] is only valid for > 1 and 
M > 0. In the case that the neutrino multiplet with a 
single species, it is necessary to set M — 0. 

The rate of accidentals should be matched to at most, 
what is feasible to follow-up with robotic optical tele- 
scopes, about 50/year. The resulting angular uncertainty 
of the multiband or multiflavor neutrino multiplet must 
also match the field of view of robotic optical telescopes, 
about 2° X 2°. Using eq. [3] we find that optical follow- 
ups are possible for 2 or more >100 GeV (which 
has already been proposed), but also (with KMSNet) 
one >100 GeV V/^i and one >100 GeV cascade. Two 
>100 GeV cascades also produces an accidental rate that 
might be appropriate, but its follow up would require a 
very large field of view telescope. Coincidences of one 
>100 GeV Vfj, and at least three DeepCore events also 
provides a good match to perform optical follow up. Ta- 
bles U and In] show the accidental false rate of neutrino 
multiplets for i/^, DeepCore events and cascades. 



Following the RMW/AB model, we have calculated the 
event expectation in DeepCore and for cascades in km^ 
detectors like IceCube and KMSNet. We find that for a 
reference core-collapse supernova at 10 Mpc ~ 4 events 
would be detected by DeepCore and ~ 6 cascades would 
be detected by km'^ detectors. These signals are strong 
enough to allow for the search of neutrinos in coincidence 
with a known supernova in the scale of 10 Mpc. In the 
case of neutrino only searches, the atmospheric neutrino 
background is higher than what we described in the text, 
as the time of the explosion can only be established with 
about 1 day precision [27| . 

A better alternative is to expand the already running 
optical follow-up programs. A single TeV z/^ event can 
provide accurate location in the sky, while a coincident 
observation with at least three 10 GeV events in Deep- 
Core would result in an acceptable false accidental rate. 
One muon event and one cascade event at >100 GeV also 
has a very low accidental rate in KMSNet, because of its 
good angular resolution for cascades. A DeepCore-like 
component in KMSNet may have better angular resolu- 
tion than in IceCube, but the level of improvement is 
limited by the kinematical direction difference between 
the neutrino and the outgoing muon or shower. 

The >10 GeV observations provide another advan- 
tage. Because they are sensitive to the pion contribu- 
tion and the >100 GeV neutrinos are sensitive to the 
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TABLE I: Expected coincident yearly rates of N events and 
M DeepCore events. We assume the z/^ events coincident in a 
circle of 1° radius and 30° for DeepCore events. We assume 
coincidence time window of 100 s. 
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TABLE IL Expected coincident yearly rates of N i/^ events 
and M cascade events. We assume the v^i, events coincident 
in a circle of 1° radius and 5° for cascade events. We assume 
coincidence time window of 100 s. 
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kaon contribution, the multiband channel helps the Ice- 
Cube/DeepCore combination to maintain sensitivity if 
actual supernovae deviate from the reference model. 

Observations of >10 GeV and >100 GeV neutrinos in 
coincidence with core collapse supernovae vi^ould be very 
strong evidence for the existence of choked jets. This 
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